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We study the superconducting transition temperature T c of the bilayer d-p model with d x 2_ y 2- 
wavelike attractive interaction based on the formalism first employed by Nozieres and Schmitt- 
Rink. In the strong coupling regime, T c obtained through this formalism are much suppressed, 
compared with those through Thouless criterion only. We also find that, whether the interlayer 
coupling exsists or not, T c is almost propotional to the Fermi energy Ef in the strong coupling 
regime. Thus, we can reproduce the essential nature in the underdoped region. 
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The pseudogap state of underdoped cuprates has been 
a very important issue in the study of high-T c supercon- 
ductors (HTSC). Since the discovery of HTSC, the pres- 
ence of an excitation gap, both in the charge and spin 
dynamics of underdoped cuprates in the normal state, 
has been indicated in several different experiments. The 
normal state transport properties (in-plane resistivity, 
Hall effect) and static susceptibility change their tem- 
perature dependences below a characteristic temperature 
T* .EP The electronic specific heat analysis has shown that 
the electronic entropy starts to decrease at the temper- 
ature well above T c u The nuclear magnetic resonance 
(NMR) experiments have found the anomalous temper- 
ature dependences both of Knight|-^bift and relaxation 
rate in almost all kinds of HTSC.Btf The suppression 
of the in- plane scattering rate well above T c has-been 
derived from the optical conductivity mesurementlf 

From all these various experimental results, we can 
deduce the existence of some kind of electronic bound 
states in the normal state of underdoped cuprates, whose 
characteristic energy should be related to T*. In order to 
explain this bound state formation, the theoretical works 
that identify T* with the singlet resonating- valence-bond 
(RVpB) formation temperture have been carried out ear- 
lier.© In the recent years, angle resolved photo-emission 
spectroscopy (ARPES) has revealed the electronic struc- 
ture of a pseudogap in the normal state, which has a 
highly anistropic momentum dependence, d x 2_ y 2 1 as well 
as in the superconducting state.LrtP 

Based on this result, it may be natural that the ori- 
gin of the pseudogap is in the preformed Cooper pair 
which cannot Bose-condensate due to its thermal fluc- 
tuation. The thermal fluctuation of the Cooper pair is 
neglected to determine T c in the original BCS theory, 
which is suitable for the description of the weak cou- 
pling regime where the electronic Fermi temperature Tp 
is much larger than the temperature T and the supercon- 
ducting carrier density is small compared to the degen- 

* E-mail: koike@ton.scphys.kyoto-u.ac.jp 



erated electronic one. However, in underdoped cuprates, 
which are believed to be in the strong coupling regime 
where the electronic Fermi energy Ep is much reduced 
by the strong correlation, the thermal fluctuation cannot 
be neglected over the wide temperature range above T c , 
where T may exceed the Fermi temperature Tp. 

For the above reason, we use the formalism first em- 
ployed by Nozieres and Schmitt-Rink (NSR),EP in order 
to estimate T c of the bilayer d-p model in such a strong 
coupling regime. Their formalism and continuum mod- 
els introduced in Ref. [)] have been adopted by the sev- 
eral authors in the past in order to study the ^jari-S.- 
dimensional strong coupling super conductivityJ13lil£2P 
The model Hamiltonian is as follows: 
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Here df ta . T (d^ ) and p x ^ ipu ) ar e the annihilation 
(creation) operator for d- and p^^-electron of momen- 
tum fe, spin a = {T,4} and layer r = {+,—}, respec- 
tively, fi is the chemical potential. The non-interacting 
part H± is represented by 
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We introduce the three types of transfer energies, t, t pp 
and t z . The first two are the intralayer transfer energies 
for the nearest neighbors between d- and p^^-orbitals 
and between p x - and p y -orbitals, respectively. The third 
is the interlayer transfer energy for between d-orbitals. 
In the following part of this letter, we take t as the unit 
of energy. The residual term H e i in eq. (1) represents 
the attractive interaction between intralayer d-orbitals, 
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described as follows: 
U d 
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where g^ = cosk x — cosk y is the d x 2_ y 2-wavelike form 
factor and N is the number of Cu sites. 

Then, we calculate the pairing correlation function for 
d-electrons_Diagonalizing H++H- by the unitary trans- 
formation,!^' we can obtain the band dispersion, e^ ±) 
and its weight for the c£-electron, , where a = {0, ±1} 
is the band index. They are defined by 
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Using these expressions, the intra-(inter-)layer pairing 
3la1 



correlation function, X+(_)j is represented by 
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Here 

fc± 



:j£ ± — /x, and /(w) represents the conven- 
tional Fermi distribution function. TVl is the number of 
fc-space lattice points in the first Brillouin zone (FBZ). 

NSR formalism takes into account the contribution 
of the phase shift from the scattering t-matrix to the 
thermodynamic potential, which leads to the equation 
that determines the critical temperature T C NSR as the 
Bose-Einstein condensation temperature for the density 
of pairs n — nt, where n is the total Fermion density and 
nf is the free part. First, the condition for the Cooper 
instability is given by 

o = i- x f(o,o) 
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Fig. 1. T ( ? /IF and T C NSR as functions of total electron number n 
in (a) for 4 Z = 0.004 and in (b) for t z = 0.204. In both panels 
the closed symbols denote the results for T < ? /IF and the open 
OIK'S for T C NSR . The circles, squares and diamonds correspond 
to Ed = 0.404, 0.80t and 1.204, respectively. The insets show 
the corresponding ratios of d-hole number njj to p-hole number 
n? for every Ed as functions of n. The correspondences of the 
symbols to Ed are the same. 
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and the equation which determines the chemical poten- 
tial /j, is described as 



n - rif (fi, T c ) = 2n h (p, T c ) 

= 4:T,9(Vqr), 



(9) 



qr 



where rjq± is a discrete pole of [1 — X± (<Z, and 
g (ui) represents the Bose distribution function. In order 
to subsequently define the Fermi energy, we introduce 
[i* as rif (p*,T c ) = n for a given total density n. Equa- 
tion (8) is called Thouless criterion, which results in the 
conventional BCS gap equation and gives the supercon- 
ducting transition temperature evaluated in the mean- 
field thory T C MF if we neglect rib and substitute n = ri[ 
in eq. (9). The summation of momentum and layer in 
eq. (9) is executed, except for the point which satisfies 
eq. (8), in order to remove the divergence. 
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Fig. 2. U d -dependence of T C NSR and T C MF for n = 4.80. The 
solid line and the dashed one denote T C NSR and T C MF respectively, 
and both are for 4 Z = 0.004 and Ed = 1.204. We obtain these 
results by the interpolation based on the n-dependence of T c for 
every U d . 



In the numerical calculation, we set e p = and control 
the filling by shifting /i. We define the Fermi energy E-p 
as follows: 



E F = fi 



E F , 



n { (Eg, 0) = 4.0, 



(10) 
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which means that Ep = i.e. /i* — E F when the highest 
band is empty. FBZ is divided into 64 x 64 mesh, and 
the other input parameters are set like, Ed = 0.40i, 0.80i, 
1.2Qf 3 t pp = 0.40i, U d = AM, t z = O.OOi for monolayer 
case or t z = 0.20t for bilayer case. When we numerically 
integrate the sum in eq. (9), we introduce the cutoff mo- 
mentum q^- — (5/32)7r assuming weak 3-dimensionality 
of real systems, in order to avoid the large fluctuation 
which is unique to lower-dimensional system. In prac- 
tice, first we investigate for a given \i, then calculate 
rib (£«, T c ) and determine /i* for the total density n. Thus, 
the nf(/x*) = n to T C NSR curve and the rif(/z) = n to T C MF 
curve can be obtained. Figure 1 shows these for the two 
diffcrnt values of t z . The difference between T C MF and 
jiNSR i ncreases mu ch near the half-filled state, which cor- 
responds ton = 5.0 for our model. While T C MF are higher 
for larger sg, T C NSR take almost the same values near the 
half-filling, independent of Ed- It means that while T C MF 
are determined by Ud/W, T C NSR are determined only by 
the filling n, where W is the bandwidth of the high- 
est band, which decreases as Ed increases. In addition, 
the suppression of T C NSR is sensitive to the exsistence of 
the interlayer transfer t z , which is supposed to introduce 
the higher dimensionality and weaken the strong thermal 
fluctuation. In Fig. 2 we compare the Ud— dependence of 
T C NSR and T C MF for n = 4.80, t z = O.OOt and e d = 1.20t. 
While T C MF increases rapidly with U d , T C NSR is almost 
saturated in the strong coupling regime. For references, 
in Fig. 3 we roughly show the momentum dependence of 
rjq±, which corresponds to the dispersion of a bosonic ex- 
citation i.e. a noncondensed Cooper pair, for two differ- 
ent fillings. rjq± becomes lower as n approaches 5.0, and 
we expect that the densities of the noncondensed Cooper 
pairs are increased and that the Bose condensation tem- 
peratures i.e. T C NSR of the pairs are much suppressed, 
compared with the temperature where the Cooper insta- 
bility occurs i.e. T C MF . In Fig. 4, we show the relation 
between T C NSR / Ep and Ep for each set of parameters. In 
our calculation, for every case T c nsr /Ef shows a max- 
imum in the intermediate coupling regime. This result 
is an artifact of approximations, and the method to re- 
move such a maximum has Jaeen already discussed by 
R. Haussmann in his paper £-f However, in the strong 
coupling regime, all qualitative results are consistent in- 
dependent of any kinds of approximations, and we find 
that T c NSR /i?F tend to settle with constant values. 

In summary, we have calculated of the bilayer 

d-p model on the assumption that d x 2_ y 2 -wavelike at- 
tractive intraction works on d-electrons. Compared with 
T C MF , T C N are much suppressed in the strong coupling 
regime. This effect is due to the thermal fluctuation of 
preformed Cooper pairs and can be weakened by the ex- 
sistence of the interlayer coupling. Near the half-filled 
state, T^ 11 is almost propotional to Fermi energy Ep. 
However, our model does not include the on-site repul- 
sive interaction among d-electrons, which is very impor- 
tant to explain the Mott-insulating state near half-filling. 
In real systems, -Ef should be much reduced near half- 
filling by this strong correlation, where T c decreases in- 
spite of strong coupling. Therefore, if we include the 
strong on-site repulsion by an appropriate method and 
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Fig. 3. q— dependence of r)q± near (0, 0) for the two different 
fillings. These data have been reconstructed from the real nu- 
merical data by an appropriate method so that their general 
view can be easily observed. The solid line and the dashed one 
denote for n = 4.1528, T C NSH = 0.01296* and for n = 4.7198, 
T NSR _ 0.12625*, respectively. Both are for tz = 0.00* and 
Ed = 1.20*. Arrows show the points for q = (q c ,0) and 

q = (o,<?c). 
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Fig. 4. The ratios of T c (= T C NSR ) to Fermi energy Ep in (a) for 
*z = 0.00* and in (b) for t z = 0.20* as functions of Ep/t. In 
both panels the open circles, squares and diamonds correspond 
to e d = 0.40*, 0.80* and 1.20*, respectively. 



renormalize Ep of the quasiparticle, we might obtain a 
more realistic picture of the underdoped cuprates. Fi- 
nally, the recent experiment has clarified that the re- 
duced gap 2Ar)/£;BT* is nearly constant for the under- 
doped Bi2Sr 2 CaCu208+,5, where 2A is the magnitude 
of the superconducting gap estimated by STS and T* is 
the temperature at which the in-plana. resistivity starts 
to deviate from the T-linear behavior £3 This result sug- 
gests the close relationship between T* and preformed 
incoherent pairs. Consequently, our future problem is 
to obtain 2Ao consistently and study the correlation be- 
tween 2A and T*(~ T C MF ) or T c (~ T C NSR ) by the use of 
more realistic interaction among (i-electrons in order to 
explain such experimental results. 
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